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ABSTRACT 

This paper presents a model for the outcome of coUisions between planetesimals in 
' a debris disk and assesses the impact of coUisional processes on the structure and 

size distribution of the disk. The model is presented by its application to Fomalhaut's 
collisionally replenished dust disk; a recent 450 /im image of this disk shows a clump 
J> ' embedded within it with a flux ~ 5 per cent of the total. The following conclusions 

I are drawn: (i) SED modelling is consistent with Fomalhaut's disk having a collisional 

CO . cascade size distribution extending from bodies 0.2 m in diameter (the largest that 

' contribute to the 850 fj,m flux) down to 7 /im-sized dust (smaller grains are blown out 

, of the system by radiation pressure), (ii) Collisional lifetime arguments imply that the 

collisional cascade starts with planetesimals 1.5-4 km in diameter, and so has a mass of 
20-30 A/q. Any larger bodies must be predominantly primordial, (iii) Constraints on 
the timescale for the ignition of the collisional cascade from planet formation models 
are consistent with these primordial planetesimals having the same distribution as the 
cascade extending up to 1000 km, resulting in a disk mass of 5-10 times the minimum 
Q . mass solar nebula, (iv) The debris disk is expected to be intrinsically clumpy, since 

$—( ' planetesimal collisions result in dust clumps that can last up to 700 orbital periods. 

, The intrinsic dumpiness of Fomalhaut's disk is below current detection limits, but 

d ■ could be detectable by future observatories such as the ALMA, and could provide the 

only way of determining this primordial planetesimal population. Also, we note that 
such intrinsic dumpiness in an exozodiacal cloud-like disk could present a confusion 
limit when trying to detect terrestrial planets, (v) The observed clump could have 
JH I originated in a collision between two runaway planetesimals, both larger than 1400 

km diameter. It appears unlikely that we should witness such an event unless both 
the formation of these runaways and the ignition of the collisional cascade occurred 
relatively recently (within the last ~ 10 Myr), however this is a topic which would 
benefit from further exploration using planet formation and collisional models, (vi) 
Another explanation for Fomalhaut's clump is that 5 per cent of the planetesimals 
in the ring were trapped in 1:2 resonance with a planet orbiting at 80 AU when it 
migrated out due to the clearing of a residual planetesimal disk. The motion on the 
sky of such a clump would be 0.2 arcsec/year, and it would be more prominent at 
shorter wavelengths. 

Key words: circumstellar matter - stars: individual: Fomalhaut - stars: planetary 
systems: formation. 
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1 INTRODUCTION 

At least 15 per cent of main sequence stars exhibit detectable 
far-IR emission in excess of that expected from the photo- 
sphere (Plets & Vynckier 1999; Lagrange, Backman & Arty- 
mowicz 2001). The spectral energy distribution (SED) of 
this excess implies that it is thermal emission from dust in 
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regions analogous to the Kuiper belt (i.e, at > 30 AU from 
the star). The few stars for which this emission has been 
imaged confirm this location for the dust and show that it is 
confined to narrow ring-Hke disks (e.g., Holland et al. 1998; 
Greaves et al. 1998; Jayawardhana et al. 1998). The short 
lifetime inferred for this dust, due to both mutual collisions 
and radiation forces, implies that it must be continually re- 
plenished, and it is thought that the dust disk is fed by the 
collisional grinding down of a population of larger planetes- 
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imals which may have formed during the system's planetary 
formation phase (Backman & Paresce 1993). Furthermore it 
has often been speculated that the absence of dust close to 
the star in these systems has been caused by clearing of the 
residual planetesimals by a planetary system. 

The debris disk images also show that the dust rings 
are neither smooth nor symmetrical: recent sub-mm images 
of the Fomalhaut disk show a clump containing 5 per cent 
of the total disk flux (Holland et al. 2002) ; there are several 
clumps in the e Eridani disk, the largest of which contains 
6-7 per cent of the disk flux (Greaves et al. 1998); a 5 per 
cent brightness asymmetry is seen in the structure of the HR 
4796 disk (Telesco et al. 2000), which could be accounted for 
by a clump containing ~ 0.5 per cent of the disk flux; Vega's 
dust disk is dominated by emission from two dust clumps 
(Holland et al. 1998; Koerner, Sargent & Ostroff 2001); and 
P Pictoris exhibits several small clumps in the mid-plane of 
the NE extension of its optical disk (Kalas et al. 2000), as 
well as a sub-mm emission peak lying ~ 34 arcsec SW of the 
star (Holland et al. 1998; Zuckerman 2001). 

While clumps are now known to be a common feature of 
debris disks, their origin, and whether this is the same in all 
cases, is still unclear. One way to form clumps is the gravi- 
tational perturbations of unseen planets, since we know that 
their effect on the orbits of planetesimals in these disks could 
cause the observed dust rings to contain both offset and warp 
asymmetries (Wyatt et al. 1999; hereafter WDT99) as well 
as resonant clumps (Wyatt 1999; Ozernoy et al. 2000). In- 
deed the observed asymmetries have been modelled by var- 
ious authors and shown to be possible evidence of planets 
of mass between IOMq and 2M,j orbiting close to the inner 
edge of the rings (typically ~ 30 AU). However, it is diffi- 
cult to test this hypothesis as the proposed planets are not 
detectable with current techniques, and such interpretations 
must be treated with caution until alternative explanations 
have been explored. 

Another mechanism for forming clumps in these disks 
is stochastic collisions between their largest planetesimals. 
A coUisionally produced disk is inherently clumpy because 
coUisional debris follows the orbit of the parent planetesi- 
mal until it has had a chance to precess around that orbit. 
The question is, do collisions of a magnitude large enough to 
produce an observable dust clump happen often enough, and 
the resulting clumps last long enough, to make a disk ap- 
pear clumpy? The possibility of observing collisional clumps 
was flrst discussed by Stern (1996) who modelled the ex- 
pected dumpiness of our own Kuiper belt, but the applica- 
tion to extrasolar debris disks has received little attention, 
as it was thought that collisions large enough to produce 
observable dust clumps would occur too infrequently, even 
if large enough planetesimals do exist in the disks. However, 
recent modelling of the asteroid belt stressed the importance 
of stochastic collisions by showing that the collisional de- 
struction of asteroids large enough to more than double the 
brightness of the asteroid belt should occur about once every 
20 Myr (Durda & Dermott 1997; see flg. 19 of Grogan, Der- 
mott & Durda 2001). This possibility may be supported by 
evidence of stochastic increases in the interplanetary dust 
flux from '^He in deep-sea sediments (Farley 1995). Also, 
dust clumps caused by collisions in the asteroid belt may 
already have been observed — the temporary dust tail ac- 
tivity of the asteroid-comet Elst/Pizarro may have resulted 



from collisions with other main belt asteroids (Toth 2000). 
This prompted Wyatt et al. (2000) to revisit the collisional 
hypothesis by applying a simple collisional model to the re- 
cently obtained image of the Fomalhaut disk (Holland et 
al. 2002). This implied that collisions were unlikely to be 
causing Fomalhaut's clump, but they concluded that this 
possibility could not be ruled out because of uncertainties 
in the disk's size distribution. 

It is the purpose of this paper to explore in depth the 
role of collisions in shaping the structure of the Fomalhaut 
disk. In particular the aim is to determine whether the clump 
observed by Holland et al. (2002) could be collisional in ori- 
gin. One reason why it is crucial to test the collisional hy- 
pothesis is that the only observation that has been proposed 
to test whether the dust clumps are associated with plan- 
etary resonances relies on the fact that such clumps would 
orbit the star with the planet — motion that should be de- 
tectable on timescales of a few years (Ozernoy et al. 2000). 
However, it may not be possible to use detection of such 
motion as an unambiguous test for the presence of planets, 
since clumps that are collisional in origin would also orbit 
the star. After summarizing the observations of the Foma- 
lhaut disk (section 2) we consider the size distribution ex- 
pected from theoretical arguments (section 3), then estimate 
this distribution observationally by modelling the observed 
SED (section 4). A model of collisions is then developed 
which is used to explore the role of collisional destruction 
in forming the observed size distribution (section 5). This 
model is then expanded to determine what magnitude of 
collision we expect to see in Fomalhaut's disk due to the 
break-up of its planetesimals (i.e., its intrinsic dumpiness; 
section 6). Finally we explore the possible scenarios in which 
collisions might have produced Fomalhaut's dust clump and 
propose another possible origin for the clump (section 7). 
The conclusions are given in section 8. 



2 THE FOMALHAUT DISK 

Fomalhaut was one of the first four main sequence stars 
found to exhibit infrared emission in excess of that from the 
photosphere and one of only a handful that have had that 
emission resolved. It is an A3V star at a distance of = 7.7 
pc that is estimated to have an age of tsys ~ 200 ± 100 Myr 
based on both evolutionary isochrone models (Lachaume et 
al. 1999; Song et al. 2001) and on the ages of the Castor 
moving group with which it is associated (Barrado Y Navas- 
cues 1998) and that of its common proper motion compan- 
ion GL879 (Barrado Y Navascues et al. 1997). The stel- 
lar parameters that have been adopted in this work are: 
Toft = 9060 K, log 3 = 4.29, M = 2Mq (Song et al. 2001); 
and L — 13Lq (Backman & Paresce 1993). In all modelling, 
the stellar spectrum is assumed to be that of a Kurucz model 
atmosphere with solar metallicity and the above parameters. 
No gas disk has been observed around Fomalhaut (Liseau 
1999). 

The photometric observations of emission from the Fo- 
malhaut disk are summarized in Table ^, and this SED is 
plotted in Fig. 0. Recently the structure of the disk was 
mapped at both 450 and 850/im using SCUBA at the JCMT 
(Holland et al. 2002; Fig. |[i). The double-lobed emission 
feature seen in the images implies that the disk is being ob- 



Collisional Processes in Extrasolar Planetesimal Disks 3 



Table 1. Photometric observations of the excess emission from 
Fomalhaut, F^, and the la errors, d_F^. The photospheric contri- 
bution has been subtracted from the observed emission using a 
Kurucz model atmosphere scaled to K = 0.992 (Bouchet, Man- 
froid & Schmider 1991). The IRAS fluxes were taken from the 
Faint Source Catalog and a colour correction applied based on 
the temperature of the excess at each wavelength. This temper- 
ature is not well known at 25 fira and the resulting uncertainty 
in the derived flux is included in the errors. The ISO fluxes were 
taken from Walker et al. (in preparation). 



A, fim 


F^, Jy 


dF,, Jy 


Reference 


25 


1.24 


0.6 


IRAS 


60 


9.80 


0.62 


IRAS 


60 


7.51 


0.92 


ISO 


80 


8.86 


0.87 


ISO 


100 


10.95 


0.66 


IRAS 


100 


11.04 


2.0 


ISO 


120 


9.77 


0.98 


ISO 


150 


7.12 


0.67 


ISO 


170 


8.18 


0.64 


ISO 


200 


3.55 


1.15 


ISO 


450 


0.595 


0.035 


Holland et al. (2002) 


850 


0.097 


0.005 


Holland et al. (2002) 



served close to edge-on and that the dust is constrained to a 
narrow ring ~ 150 AU from the star. This is consistent with 
the disk's SED, which implies emission from cool dust at a 
single temperature, as well as with previous observations of 
the disk's structure (Holland et al. 1998; Harvey & JefTrys 
2000). The images at both wavelengths also show an asym- 
metry in the brightness distribution. The spatial structure of 
the disk seen in the 450 ^m image (Fig. ^) was modelled in 
Holland et al. (2002), where they showed that the observa- 
tion is consistent with emission from a smooth axisymmetric 
ring (Fig. ^) embedded within which is a bright clump with 
a flux equal to ~ 5 per cent of the total disk flux (~ 30 mjy; 
Fig. ![:). Holland et al. (2002) argue that this clump is real 
and associated with the star, as the asymmetry is signifl- 
cantly above the noise at both wavelengths, and there is a 
low probability of contamination from a background source. 

In the modelling presented here, we use the model pre- 
sented in Holland et al. (2002) for the smooth ring. This 
has a radial distribution of cross-sectional area a{r) - where 
a{r)dr is the fraction of the total cross-sectional area crtot 
between r and r+dr - that is indicative of a ring that has 
a mean radial distance of ~ 150 AU and a width of ^ 50 
AU. This model is however described by the distribution of 
the orbital elements of the material in the ring (see, e.g., 
WDT99). Material in their model has average orbital incli- 
nations and eccentricities of / = 5° and e — 0.065. While 
this modelling could not uniquely constrain these parame- 
ters, these values result in a best fit to the data and are 
consistent with indications later in this paper that we are 
observing a collisional cascade for which such high values 
are likely (section 5.3). 



THEORETICAL PLANETESIMAL SIZE 
DISTRIBUTION 
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Figure 1. The SED of emission from Fomalhaut after subtraction 
of the stellar photosphere (see Table |l|). The asterisks, diamonds 
and triangles correspond to the IRAS, ISO and SCUBA fluxes, 
respectively. The lines represent model flts to these data (see sec- 
tion 4) that assume that the disk's spatial distribution is deflned 
by a fit to the 450 fim image (Holland et al. 2002), and that it is 
comprised of solid grains of 1/3 (by volume) amorphous silicate 
and 2/3 organic refractory material (e.g., Li & Greenberg 1997, 
1998). The grains' size distributions in the models extend down to 
the radiation pressure blow-out limit, Dmin = 7 Mui, with power 
law indices given by q^. 



3.1 Basic idea 

Debris disks are thought to have formed in the first 10 Myr 
or so of a star's life through canonical planet formation 
processes, wherein collisions between initially sub-^m dust 
particles resulted in their coagulation and eventual growth 
into planetesimals and maybe even planets (e.g., Lissauer 
1993). Toward the end of this process the collisional process- 
ing was reversed whereupon collisions between planetesimals 
more frequently resulted in their net destruction rather than 
growth. Such a change is likely to have been the consequence 
of an increase in the average relative velocity of collisions, 
probably caused by some combination of the removal of the 
gas disk (which was acting to damp the eccentricities and in- 
clinations of disk material) and gravitational perturbations 
from planets or large planetesimals that formed in the disk 
(Kenyon & Bromley 2001), although this cross-over stage 
is not well understood. In the debris disk systems, we know 
that by the end of this process, regardless of whether planets 
or planetesimals formed close to the star, a disk of planetes- 
imals was able to form and subsequently survive in a ring 
at > 30 AU from the star (Lagrange et al. 2001). Once their 
collisional velocities are pumped up, these planetesimals col- 
lide and break up into smaller fragments, initiating a colli- 
sional cascade. A typical planetesimal in such a cascade is 
most likely to have been created by the break-up of a larger 
parent body in a collision with another large planetesimal. 
This parent body would most likely have been created by the 
break-up of an even larger body, and the planetesimal itself 
will most likely end up as a parent body for planetesimals 
smaller than itself. It is the large quantities of small dust in 
the cascade, and its commensurate large surface area, which 
makes a debris disk observable. This dust is replenished as 
long as there remains a supply of large primordial bodies to 
feed the cascade. 
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Figure 2. 450 fim images of tlie Fomalliaut disk (adapted from 
Holland et al. 2002): (a) observation taken with SCUBA at the 
JCMT, (b) axisymmetric smooth disk model, and (c) the residuals 
(= observation - model), which show that the asymmetry in the 
observation could be explained by a bright clump embedded in 
the smooth disk. All contours are spaced at la = 13 mjy beam"'^ 
levels. For (a) and (b) the lowest (solid white) contour is set at 
3fT, while for (c) the solid white line indicates the zero level and 
the black solid and dashed lines are the positive and negative Icr 
levels respectively. The dashed white oval shows the approximate 
inner edge of the iriid-plaiie of the disk, a 125 AU radius ring 
inclined at 20° to the line of sight. The coordinates are given as 
their offsets from the position of the star. The stellar photosphere 
has been subtracted from the observation. 



3.2 Size distribution of an infinite collisional 
cascade 

The equilibrium size distribution resulting from a collisional 
cascade can be found from theoretical arguments (Dohnanyi 
1969): 

n(D)ocD^-^^^ (1) 

where n{D)d.D is the number of planetesimals of size be- 
tween D and D+dD. For a self-similar collisional cascade it 
can be shown that ga = 1.833 (Tanaka, Inaba & Nakazawa 
1996), however if the strength of a planetesimal is size de- 
pendent (see section 5), then the slope of the equilibrium 
distribution could be slightly shifted from 1.833 (see fig. 12 
of Durda & Dermott 1997; hereafter DD97). 

3.3 Size distribution of small grains 

The collisional cascade distribution is only expected to hold 
for disk particles that are large enough not to be affected 
by radiation forces (WDT99). Particles smaller than a few 
/im are blown out of the system by radiation pressure on or- 
bital timeseales and all particles spiral into the star due to 
Poynting-Robertson (P-R) drag on timeseales that are pro- 
portional to their size. It can be shown, however, that the 
high densities of the observed extrasolar disks result in colli- 
sional lifetimes for dust of all sizes that are shorter than their 
P-R drag lifetimes and so the latter effect can be ignored 
(WDT99). Thus here we make the assumption that the size 
distribution holds down to the radiation pressure blow-out 
limit Dmin (i.e., those for which /3 = Frad/^grav > 0.5), be- 
low which the distribution is cut-off. 

A sharp cut-off below -Dmin is, however, a simplification. 
For a start, a sharp cut-off would cause a "wave" in the equi- 
librium size distribution (Durda, Greenberg & Jedicke 1998; 
hereafter DGJ98), since in an equilibrium situation particles 
close to the blow-out limit should have been broken up by 
the missing population of particles smaller than themselves. 
Since these are not removed their number increases and this 
affects the particles they are expected to break up, and so 
on. Also, while it may be possible to ignore the contribu- 
tion of particles smaller than the blow-out limit to the total 
cross-sectional area of disk material (e.g., if their lifetime is 
much shorter than their collisional lifetime, WDT99), the 
size distribution of gravitationally bound grains close to the 
blow-out limit could be affected by their destruction in col- 
lisions with the blow-out grains (Krivov, Mann & Krivova 
2000). The size distribution of these bound grains would 
be further complicated by radiation pressure which causes 
them to have large eccentricities, thus reducing their num- 
ber density in the disk region (see, e.g., fig. 4 of Augereau 
et al. 2001). 

3.4 Size distribution of large planetesimals 

At the large size end, planetesimals that have a collisional 
lifetime that is longer than the age of the system do not form 
part of the collisional cascade, and their distribution must 
be primordial. However it is unclear what this would be, as 
even the primordial distribution of solar system objects is 
unknown. It is uncertain whether the distribution of aster- 
oids larger than 30 km is primordial due to uncertainties 
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in the collisional properties of these asteroids (DGJ98). The 
distribution of Kuiper belt objects larger than 100 km, which 
appear to have a relatively steep distribution {qa = 1.83- 
2.27, Luu & Jewitt 1998; Gladman et al. 1998), is also dif- 
ficult to interpret, since while the objects themselves might 
be primordial, their distribution may not be. The sweep- 
ing of Neptune's resonances through the young Kuiper belt 
as Neptune's orbit expanded early in the history of the so- 
lar system (Hahn & Malhotra 1999) would have resulted in 
a size distribution of objects trapped in resonance that is 
steeper than the primordial distribution (Jewitt 1999). The 
primordial distribution expected from theoretical grounds 
is also unclear, since simulations of planet formation result 
in distributions that are usually steeper than the cascade 
distribution, but with anywhere between 2 and 2.5 (e.g., 
Kokubo & Ida 2000; Kenyon 2002). In any case, the size 
distribution of the largest planetesimals in the distribution 
is likely to be dominated by the runaway growth process, 
whereby these grew faster than the rest of the population, 
causing a bimodal size distribution (e.g., Wetherill & Stew- 
art 1993; Kokubo & Ida 1996). Runaway growth is discussed 
further in section 7.2. 



by their porosity, p, and the volume fraction of the gap in 
the grains that is filled with ice, gicc. We use the optical con- 
stants for each of the components taken from Li & Greenberg 
(1997, 1998) and compute those of the composite material 
using Maxwell-Garnett Effective Medium theory (Bohren & 
Huffman 1983). The absorption efficiencies of these grains, 
Qabs(A, D), are then calculated using Mie theory (Bohren & 
Huffman 1983), Rayleigh-Gans theory or geometric optics 
in the appropriate limits (see Laor & Draine 1993). These 
then define the temperature of the grains, T[D, r), and their 
radiation pressure blow-out limit, Dmin (see, e.g., WDT99). 

Thus the SED model is defined by the following four 
parameters: atot (the total cross-sectional area in the disk), 
gd, p, and gicc. The flux in Jy from the disk at a given wave- 
length is given by (e.g., WDT99) 



2.35 X 10""ii^Vto 



Qab.(A,D)a(D) 



B^[T{D,r)]G{r)ArAD, 



(2) 



where (j{D) = 0. 25tv n{D) /atot and utot is in AU^. Here 
we have assumed that the disk is optically thin at all wave- 
lengths and that there is no disk self-heating. 



4 SED MODELLING 
4.1 The SED model 

A disk's SED is determined by three factors: the spatial 
distribution of material in the disk, its size distribution, and 
its optical properties (determined by its composition) . While 
these factors could be interrelated, we assume here that they 
are independent. For the spatial distribution we assume that 
particles of all sizes have the spatial distribution inferred 
from the modelling of the 450 fim images defined by a-{r) 
(see section 2). For the size distribution we consider one 
that follows equation (|l|), and assume that this distribution 
is cut-off below Dmin, the radiation pressure blow-out limit, 
and extends up to large enough sizes that this limit does not 
affect the SED. 

For the particles' composition we use the core-mantle 
model developed by Li & Greenberg (1997) for interstel- 
lar dust, which has since been applied to cometary dust 
(Greenberg 1998) and to dust grains in extrasolar disks (Li & 
Greenberg 1998; Augereau et al. 1999). It seems reasonable 
to apply this model to extrasolar disks, since these grains 
are likely to be reprocessed (albeit significantly) interstellar 
grains. In this model the dust is assumed to be aggregates 
of core-mantle grains with a silicate core (p = 3500 kg m~"^) 
and an organic refractory mantle comprised of UV photo- 
processed ices which accreted onto the silicate cores in the 
interstellar medium (p = 1800 kg m~^). These aggregates 
are assumed to be porous and to have a fraction of wa- 
ter ice (p = 1200 kg m~'^) filling the gaps; this is ice that 
would have frozen onto the growing grains in the protoplane- 
tary disk. Since interstellar grains are amorphous and would 
only be crystallized at the high temperatures expected close 
to stars, we assume that both the silicate and ice compo- 
nents of the grains are amorphous. We also assume a vol- 
ume fraction of silicates to organic refractory material in the 
core-mantle of 1:2 as found for interstellar-type grains (Li & 
Greenberg 1997). Thus the particles' composition is defined 



4.2 The modeUing procedure 

We determined the best fit for each particle composition 
(defined by p and gico) by minimizing = [('P'i'.,obs - 

Fvi. mod) /dF„^ ^ohs]^ , where atot was set to scale the SED to fit 
the 850 pm flux. Fig. |^ shows how this procedure was applied 
for solid grains (for which p = 0, Dmin = 7 pm, p = 2370 kg 
m"'', and which by definition contain no water ice from the 
protoplanetary disk). The best fit gave q^ — 1.84, atot ~ 33.7 
AU^, and = 43.4. This type of modelling constrains gd 
very well since it has a large effect on the SED: increasing 
5d increases the contribution of small grains and so results 
in an increase in the flux shortward of about 60 pm (due to 
the higher temperatures of small grains) but a decrease in 
the fiux at longer wavelengths (since small grains emit less 
efficiently than larger ones at long wavelengths)^ 

Increasing gd also steepens the sub-mm spectral slope, 
defined by a, where oc A~", because the emission here 
is dominated by grains of a size comparable with the wave- 
length so that the ratio of fluxes at two wavelengths is de- 
termined by the ratio of cross-sectional area in grains of size 
comparable to those wavelengths. This effect can be quan- 
tified by a simple model in which all the dust is at the same 
distance from the star and has the same temperature, T, 
and has emission efficiencies that are given by 



Oabs(A,D) 



1 for A < D 

(X/D)"- for A > D 



(3) 



Putting this into equation (P), we find that in the Rayleigh- 
Jeans regime: 



F^ <x 



X^-Sqa foj. n > 3gd - 5 
A-2-" for n < 3gd - 5; 



(4) 



^ It is not immediately obvious from Fig. |l| that this cross-over 
occurs at 60 pm because, as these models have been scaled to 
-Fi^(850 pm), those with a higher gd S'lso have a higher ctot. 
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Figure 3. The best fits to the observed SED for five models with 
diff'crcnt grain compositions (see legend). The resulting best fits 
for the size distribution, the total cross-sectional area in the disk, 
and the for the five models (top to bottom in the legend) arc: 
ijd = [1.84,1.81,1.87,1.81,1.88], crtot = [33.7,33.2,27.1,34.7,23.4] 
AU^, and = [43.4,75.3,101,118,117]. The minimum grain size 
in these models is: D^\^ = [7,13,9,124,11] /xm. 



i.e., increasing steepens the slope up to a constant value 
for > 5/3 + n. Thus the sub-mm spectral slope is largely 
determined by the size distribution of particles of sub-mm 
sizes. This demonstrates the importance of an SED model 
not only giving a good fit to the whole SED, but also 
giving a good fit to Oobs = 2.85. In this simple model, the 
observed slope implies that either — 1.95 and n > 0.85, or 
n = 0.85 and > 1.95. In fact we find that for solid grains 
ce — aoba ~ 4.9(gd — 1.86) and so qd — 1.84 fits the sub- 
mm slope relatively well. The discrepancy from the simple 
model is because real grains have emission features that are 
not reproduced by equation rtS). 



4.3 Modelling results 

The analysis in Fig. |^ was repeated using grains with a range 
of porosities and ice fractions, but found that increasing ei- 
ther p or gico resulted in a worse fit to the data (i.e., a higher 
X^, see Fig. ^). The problem with increasing the porosity of 
the grains is that even though solid grains are less affected 
by radiation pressure than are porous grains (mostly due to 
their higher density), and so can exist in the disk down to 
much smaller sizes (see Dmin in the caption of Fig. the 
disk's smallest porous grains are still hotter than its smallest 
solid grains because they appear to be made up of an ag- 
glomeration of much smaller particles. Thus the SED result- 
ing from a porous-grained disk peaks at shorter wavelengths 
than a solid-grained one (see Fig. Higher porosity grains 
also have fewer emission features which results in a flatter 
sub-mm slope. The problem with increasing the ice fraction 
in the grains is that they have a range of emission features 
in the far-IR, but not of a shape that increases the sub-mm 
slope to the same levels as for non-icy grains. The resulting 
spectrum is too peaked in the 40-90 fim wavelength region 
(see Fig. and because icy grains are inefficient absorbers 
of starlight, their cool temperature results in a poor fit to 
the mid-IR flux. A poor fit to the 25 ^m flux would not in 
itself rule out a model, since there may be a contribution 
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Figure 4. The temperatures of different sized dust grains at dif- 
ferent distances from Fomalhaut, assuming these are solid grains 
composed of 1/3 (by volume) amorphous silicate and 2/3 organic 
refractory material. 



to the observed flux from an, as yet undetected, hotter disk 
component such as that inferred to exist in the HR 4796 
system ( Augereau et al. 1999) . We conclude that while solid 
non-icy grains are the most likely candidates for the grains 
in this disk, SED modelling cannot set strict constraints on 
the composition of the grains without a more accurate de- 
termination of the disk's mid-IR emission. However, a more 
robust feature of the SED modelling is the size distribution 
parameter q^, which lies in the range 1.81-1.88 for all mod- 
els. 

In the rest of this paper we use the parameters (fmin, p, 
gd and atot ) derived from this modelling assuming solid non- 
icy grains, since this composition is also consistent with that 
expected on theoretical grounds. Grains that grow through 
grain-grain collisions in a protoplanetary disk are expected 
to have significant porosity (e.g., Wurm & Blum 1998), simi- 
lar to that of grains resulting from the sublimation of comets 
in the solar system (p — 0.85-0.95). Grain-grain collisions, if 
energetic enough however, can lead to significant grain com- 
paction (Dominik & Tielens 1997). Thus we expect grains 
produced by comet sublimation to have a high, primordial- 
like, porosity (such as has been inferred for the grains around 
around f3 Pictoris and HR 4796A; Li & Greenberg 1998, 
Augereau et al. 1999), whereas grains at the end of a colli- 
sional cascade should be more compact, maybe with a poros- 
ity close to that of stone meteorites {p — 0.05-0.3; Flynn, 
Moore & Klock 1999). The size distribution we have inferred, 
qd close to 1.84, is consistent with the dust having its origin 
in a collisional cascade (section 3.2); collisional lifetime argu- 
ments in the following section also support this hypothesis. 
Also, we expect little contribution to the dust population 
from cometary sublimation because of the low temperatures 
in the disk (see Fig. Thus a low level of porosity is to be 
anticipated for the dust. 

It is harder to predict the amount of water ice on the 
grains. Certainly the low temperatures in the disk imply that 
while the grains were forming (i.e., before they were com- 
pacted by collisions), water vapour should have condensed 
onto the aggregates, at least partially filling the gaps in the 
grains (Pollack et al. 1994; Li & Greenberg 1998). However, 
this ice may have since been sublimated by repeated col- 
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Figure 5. The contribution of different sized dust grains to the 
observed fluxes for the best fit model to the SED (i.e., that shown 
in Fig. ^with = 1.84). The area under the curves between two 
grain sizes defines the proportion of the flux that originates from 
this size range; the total area across all sizes is unity. 



lisional heating, or removed by sputtering by stellar wind 
ions (Jurac, Johnson & Donn 1998). Or indeed this ice may 
still exist on the grains, yet it could look more like the or- 
ganic refractory material in our model due to the formation 
of an irradiation mantle by stellar UV photons and galactic 
cosmic rays (e.g., Johnson et al. 1987; Cooper, Christian & 
Johnson 1998; Jewitt 1999). We note that water ice features 
have been detected in some, but not all, Kuiper belt objects 
(e.g., Luu, Jewitt & Trujillo 2000; Brown, Blake & Kessler 
2000). 

4.4 Limitations for large planetesimals 

The size distribution inferred from the SED modelling can 
only be applied to grains that contribute to the observed 
SED fluxes. Fig. |^ shows that planetesimals larger than ~ 
0.2 m in diameter contribute less than 5 per cent of the 850 
yitm flux (and even less at shorter wavelengths). Thus we can 
only derive the size distribution for grains that are smaller 
than ~ 0.2 m (although, as we will see in the next section, 
we can extend this up to larger sizes by collisional lifetime 
arguments). Since the total mass of the Fomalhaut disk, 
assuming that the size distribution with qd = 1.84 extends 
up to Dmax, is givcu by 

mtot/Mffi ^ 0.5DS;S, (5) 

where Dmax is in m, this means that we see a dust mass of 
~ 0.2M®. 

5 CATASTROPHIC COLLISIONS 

The outcome for a planetesimal of size D when it is impacted 
by another of size Dim is determined by the specific incident 
kinetic energy, 

Q = 0.5{Din./Dfvli. (6) 

The shattering threshold, Qg, is defined as the specific inci- 
dent energy required to break up the planetesimal, resulting 
in fragments with a range of sizes up to half the mass of 



the original planetesimal. Collisions with Q < Qg result in 
cratering whereby some material is ejected, but the plan- 
etesimal remains largely intact. More energetic collisions re- 
sult in the planetesimal being reduced to more and smaller 
fragments^ In general larger planetesimals have lower shat- 
tering thresholds, since they contain larger flaws which are 
activated at lower tensile stresses (Housen & Holsapple 1990, 
1999). Very large planetesimals, however, are strengthened 
by their gravitational self-compression which inhibits crack 
propagation (Davis et al. 1985). After the collision, the frag- 
ments disperse due to the kinetic energy imparted to them 
in the collision. For the break-up of a large planetesimal 
(D > 150 m), however, this energy may not be enough to 
overcome the gravitational binding energy and some frag- 
ments may reaccumulate into a rubble pile (e.g., Campo 
Bagatin, Petit & Farinella 2001; Michel et al. 2001). A catas- 
trophic collision is defined as one in which Q > Q^, where 
the dispersal threshold, Q^, refers to an impact in which 
the largest fragment resulting from the collision after reac- 
cumulation heis taken place (i.e., this fragment could be a 
rubble pile) has half the mass of the original planetesimal. 
The break-up of a small planetesimal (i.e., D < 150 m), for 
which ~ Qg, is said to occur in the strength regime, 
while the break-up of a larger planetesimal is said to occur 
in the gravity regime, where 3> Qs- 

A planetesimal's catastrophic collisional lifetime de- 
pends on the disk's size distribution, since only collisions 
in which Dim > Dtc = AtcD have enough energy to be 
destructive, where 

Ate = {2Q*u/viy^\ (7) 

We use the term threshold-catastrophic to denote a collision 
between planetesimals of size D and Dtc(D), and wtc(D) is 
the relative velocity of such collisions. Here we assume that 
all collisions have the same average relative velocity at large 
separations^, 

Wrcl = f(e,I)v]^, (8) 

where /(e, /) is a function of the average eccentricities and 
inclinations of the planetesimals given as Vl.25e2 (Lis- 
sauer & Stewart 1993; Wetherill & Stewart 1993), and vy^ 
is the keplerian velocity at this distance from the star. In 
the Fomalhaut model /(e, 7) ~ 0.11; thus v^ci ~ 0.4 km s~^ 
at the mean distance of 150 AU. The relative velocity of a 
collision is then 

2 For collisions with very small relative velocities there are two 
other possible collisional outcomes: the planetesimals could re- 
bound in an inelastic collision, or they could stick together. While 
these possibilities are important in protostellar disks, they are not 
considered in this paper, since the collisional velocities in debris 
disks are thought to be too high for these modes to be important. 

The distributions of the eccentricities and inclinations of plan- 
etesimals predicted from planet formation models are depen- 
dent on their size. However, in a debris disk undergoing a col- 
lisional cascade these distributions would depend on the mecha- 
nism which causes the high eccentricities and inclinations in these 
disks; e.g., if this mechanism is the continual gravitational stirring 
of massive planetesimals or planets in the disk, these distributions 
would be uniform with planetesimal size, whereas if this mech- 
anism was stirring by a passing star these distributions would 
evolve with time after the event and would be size dependent 
(Kenyon & Bromley 2002). 
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viol = t'rcl + t'osc(-D, Am), (9) 

where the mutual escape velocity of the two planetesimals 
is given by 

vl^D, A„0 = {2/i)-KGp^±^. (10) 

The increase in impact velocity due to gravity (i.e., gravita- 
tional focussing) becomes important for collisions for which 
Vesc > t'roi, i.e., for planetesimals with D > v^ci/y^ (2/3)7rGp. 
For Fomalhaut this corresponds to planetesimals larger than 
about 700 km in diameter, or 0.6 per cent of a lunar mass. 

5.1 Dispersal threshold model 

Several authors have studied how Qg and Qd vary with 
planetesimal diameter for a variety of materials (e.g., 
rocks/ice) with a range of structural properties (e.g., homo- 
geneous/rubble pile), and for collisions with a range of rela- 
tive velocities (e.g., km s~^ debris disk-like/m s~^ protostel- 
lar nebula disk-like impacts) and impact parameters (e.g., 
head-on/glancing blow). The techniques that have been used 
to study collisions range from lab experiments which pro- 
vide a direct measure of a collisional outcome (e.g., Fuji- 
wara et al. 1989; Housen & Holsapple 1999), to theoreti- 
cal studies (e.g., Petit & Farinella 1993), to interpretation 
of the distributions of the asteroid families (e.g., Tanga et 
al. 1999; Cellino et al. 1999) or the main-belt population 
(DD97; DGJ98), to computational modelling using smooth 
particle hydrodynamics (SPH; Love & Ahrens 1996; Benz 
& Asphaug 1999, hereafter BA99). However, most of these 
studies are valid only in specific regimes (e.g., in the strength 
or gravity regimes). 

In this paper we use the results of BA99 who use SPH 
to model the fragmentation of both solid ice and basalt of 
a broad range of sizes (up to 200 km diameter) for impacts 
with a range of impact parameters and with relative veloc- 
ities comparable to those expected in extrasolar planetesi- 
mal disks. In the rest of this paper Qq refers to the disper- 
sal threshold averaged over all impact parameters (BA99). 
The BA99 models reproduce the well-known features of the 
Qd vs D plot; i.e., the dispersal threshold decreases with 
size (oc 75-(0-36-o.45) -J ^j^g strength regime due to the 
lower shattering strength of larger planetesimals, but in- 
creases with size (oc jji is-i-^s-j -j-j^g gravity regime due to 
the extra energy required to impart enough kinetic energy 
to the collisional fragments to overcome the planetesimal's 
gravity. However, the impact strengths they derive for cm- 
sized ice grains impacted at ~ 3 km s~^ are 20-50 times 
higher than those directly measured in the laboratory (e.g., 
Arakawa 1999; Ryan, Davis & Giblin 1999). Their resuhs 
for basalt impacted at 3 km s~^ do not show the same dis- 
crepancy (Holsapple 1994). 

While the SED modelling suggested that there is little 
water ice in the Fomalhaut grains, the structural properties 
of the organic refractory mantle would be close to that of the 
ice studied in BA99, since the mantle is essentially comprised 
of UV photoprocessed ices. Thus the collisional properties 
of planetesimals in Fomalhaut's disk would lie somewhere 
between those of rock and ice. Here we consider three dif- 
ferent models for Q^: Ice refers to the BA99 result for ice 
impacted at 0.5 and 3 km s~^; Weak Ice also refers to the 
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Figure 6. The collisional properties of planetesimals of differ- 
ent sizes in the Fomalhaut disk for the three models (Ice, Weak 
Ice, and Basalt) based on the SPH modelling of Benz & Asphaug 
(1999): (a) the dispersal threshold; (b) the relative size of an 
impactor required to catastrophically break up a planetesimal. 
The dispersal threshold decreases with size in the strength regime 
(D < 150 m) due to the lower shattering strength of larger plan- 
etesimals, but increases with size in the gravity regime due to 
the extra energy required to impart enough kinetic energy to the 
collisional fragments to overcome the planetesimal's gravity. The 
relative size required for a catastrophic impact varies with D as 
does Qd except that it flattens off for D > 700 km because the 
extra energy required to disperse the collisional fragments is pro- 
vided by the higher impact velocities, which have been enhanced 
by the planetesimals' gravitational attraction (from '^0.4 km s^^ 
at 100 km to ~ 0.8 km s"^ at 1000 km). 



BA99 results for ice, except that the part of a planetesimal's 
Qd resulting from its strength is reduced by a factor of 50 to 
reflect the discrepancy with laboratory results; and Basalt 
refers to the BA99 result for basalt impacted at 3 and 5 km 
s~^. In all cases we have used p = 2370 kg m~'^ and the 
BA99 results at the two impact velocities have been scaled 
to the relative velocity of a threshold-catastrophic collision, 
iitc, assuming that Qd oc I'Jof This procedure is iterative 
since to calculate Qd we need to know iitc, for which we need 
to know Xtc, which is itself determined by both Qd and Wtc. 
The resulting models for Qd and Xtc are shown in Fig. ^. 

While the BA99 models apply only to planetesimals up 
to 200 km in diameter, we have assumed that their results 
can also be scaled up to 1000 km diameter. We do not ex- 
pect their models to be applicable much beyond this, since 
in this regime collisions occur at the planetesimals' escape 
velocity. As some energy is lost in the collision, only a small 
fraction of the collisional fragments can escape, and so such 



Collisional Processes in Extrasolar Planetesimal Disks 9 



collisions must result in net accretion rather than destruc- 
tion (i.e., Qd — » oo). Such an effect would not have been 
picked up in BA99, since their experiments were all carried 
out with impact velocities that are more than an order of 
magnitude under the escape velocities. Collisions between 
such massive bodies, which further we expect to be rubble 
piles, are discussed in section 7. 



5.2 Catastrophic collision timescale 

The rate of impacts from planetesimals in the size range Dim 
to -Dim+ dDim falling onto a planetesimal of size D is given 
by i?coi(D, Am)dAm, where (e.g., Opik 1951) 



Planetesimal mass, M^a 



Rcol{D, Am) = f{D, Am)cr(r, e, 4>)V,,1. 



(11) 



In this expression, (T{r,9,(f)) is the average cross-sectional 
area density of disk material, which, for a ring of radius r, 
radial width dr, and that has a (small) mean orbital incli- 
nation of particles in the ring of /, is approximately 



a{r,9,(j)) = (Ttot /(47rr- dri). 



(12) 



Also, /(D, Am)f7"(r-, 6', (/))dAm is the cross-sectional area 
density of planetesimals in the given size range that the 
planetesimal of size D sees, where 



f{D, Am) = a(A. 



D + Din 
Am 



1 



3(D,An 



.(13) 



The last term in equation ( p^ is the gravitational focussing 
factor, which becomes important for collisions with plan- 
etesimals in the Fomalhaut disk that are larger than about 
700 km. 

A planetesimal's catastrophic collision timescale is thus 
given by 



tcciD) 



: tper(rdr/atot)[2///(e, /)]//cc(D), 



(14) 



where fpor is the orbital period at r, and fcc{D) is defined 
by 



/cc(D) = 



f{D, Am)dAn 



(15) 



where Ac is the larger of XtcD or Dmin. 

Fig. shows the catastrophic collisional lifetime of plan- 
etesimals of different sizes in the Fomalhaut disk for the 
three models of their collisional properties, using the size 
distribution inferred from the SED modelling (qd = 1.84) 
assuming that this is truncated at 1000 and 10,000 km. 
Normally a planetesimal is broken up by one that just has 
enough energy to do so (i.e., in a threshold-catastrophic col- 
lision), since of all the planetesimals that can destroy it, 
these have the largest total cross-sectional area in the disk. 
However, the enhanced cross-sectional area of gravitation- 
ally focussing planetesimals can mean that a planetesimal 
is more likely to end its life by being accreted onto one of 
these than by being destroyed by a smaller planetesimal, 
hence the shorter collisional lifetime of large planetesimals 
when -Dmax = 10, 000 km. As there is less cross-sectional 
area available to destroy larger planetesimals, there is a 
steady increase in collisional lifetime with size; e.g., a simple 
model which assumes that most of this area is in planetesi- 
mals of size XtcD, where Xtc < 1, shows that if Qd oc 




Figure 7. The collisional lifetime of planetesimals of different 
sizes in the Fomalhaut disk. The solid, dotted and dashed lines 
are for the three models for the collisional properties of these 
planetesimals described in the text. The two lines for each model 
assume that the size distribution inferred from the SED modelling 
(gj = 1.84) extends up to planetesimals of size 1000 and 10,000 
km (distributions with larger planetesimals have lower collisional 
lifetimes). 



then tec oc D(3+<')9d-{5-t-a)^ implying that the BA99 mod- 



els should have tec oc D 

tec oc 



in the strength regime and 
in the gravity regime, consistent with Fig. ^ 



5.3 Implications for Fomalhaut's size distribution 

While the collisional lifetime of all planetesimals contribut- 
ing to the Fomalhaut observations {D < 0.2 m) is model 
dependent, because each requires a collision with different 
sized planetesimals for catastrophic destruction (Fig. |^), 
all models predict a lifetime that is shorter than the age of 
the system (200 Myr) by more than an order of magnitude 
(Fig. mj . Thus the primordial populations of this size range of 
particles has already been substantially depleted. However, 
these populations will also have been replenished with colli- 
sional fragments from the break-up of larger planetesimals. 
Since the largest particles that we can see are replenishing 
the smaller particles, tt is an unavoidable conclusion that 
this part of the size distribution forms a collisional cascade. 
The real questions are: (i) is the cascade in equilibrium, and 
(ii) how far does it extend up in size, or in other words, what 
is feeding the cascade? 

In answer to (i), while the parameterization of the 
size distribution assumed in the SED modelling is only 
an approximation (e.g., sections 3.2 and 3.3), the derived 
slope, g(d) = 1.84, is very close to that expected for an 
equilibrium cascade. Also, the short collisional lifetimes in 
Fig. 1^ imply that the cascade should have reached equilib- 
rium within the age of the system. 

In answer to (ii), the short collisional lifetime of the 
0.2 m diameter planetesimals we see today (0.4-9 Myr de- 
pending on their composition. Fig. ^) implies that if there 
were no larger planetesimals with which to replenish their 
number, then while all 0.2 m planetesimals we see must be 
primordial, their current population should only be a small 
fraction of their original population. In fact, it must be an 
infinitely small fraction of the original population, since the 
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coUisional lifetime of these planetesimals would have been 
even shorter at earlier epochs when their population was 
more massive (this lifetime is approximately inversely pro- 
portional to its mass). In other words, unless the coUisional 
cascade was only initiated in the last few Myr, the large 
number of 0.2 m planetesimals we observe today implies 
that there must be a population of larger bodies. Similar ar- 
guments can be applied for all planetesimals for which their 
coUisional lifetime is shorter than the time since collisions 
became destructive (which we assume here is approximately 
the age of the system); the frequency of collisions amongst 
these large bodies means that their distribution would also 
be expected to follow that of a coUisional cascade. Thus we 
infer that the coUisional cascade starts with planetesimals of 
a size for which their coUisional lifetime is equal to the age 
of the system, i.e., those with -Dmax(cc) = 1-5-4 km (Fig. m), 
implying that the mass of material in the coUisional cas- 
cade is TUcc ~ 20-30 (equation |^. If we had assumed 
that collisions became destructive more recently than 200 
Myr ago, we would have inferred the cascade to start with 
smaller planetesimals. However, we note that a few km is 
also the size of planetesimal thought to have fed the coUi- 
sional cascade in the young (< 1 Gyr) Kuiper belt (Kenyon 
2002). 

We also predict that any planetesimals larger than 1.5-4 
km are predominantly primordial and so their size distribu- 
tion is also primordial; these will have contributed little to 
the coUisional cascade. Their size distribution can be some- 
what constrained from considerations of the total mass of 
the disk, which is dominated by that of the largest bodies. 
If the primordial size distribution connects smoothly to the 
coUisional cascade distribution, and can be defined by the 
power law exponent qp between -Dniax(cc) and Dmax, then 
the total mass of primordial material would be 

mp/mcc = 0.48(6 - 3gp)-'[(Dn.ax/Anax(cc))'"""' - 1]. (16) 

The minimum mass solar nebula (Weidenschilling 1977), if 
it extended out to such large distances, would have had a 
mass of solid material of 0.1-0.2 Mj in the region 125-175 
AU. Planet formation models predict that for the gravita- 
tional perturbations of large planetesimals in the disk to 
have stirred the disk sufficiently to ignite a coUisional cas- 
cade within 200 Myr, the planetesimals in Fomalhaut's disk 
would have had to have grown to at least 1000 km in di- 
ameter and the disk must be at least 10 times more mas- 
sive than the minimum mass solar nebula (see fig. 5 of 
Kenyon & Bromley 2001). Thus a reasonable distribution 
for the primordial planetesimals, and one which we will use 
for reference later in the paper, is one with gp = 1.84 and 
fmax = 1000 km, since such a disk would have a mass of 
~ 1.1 Jupiters, i.e., 5-10 times that of the minimum mass 
solar nebula. We note that 1000 km is similar in size both 
to the largest asteroids and Kuiper belt objects (Trujillo et 
al. 2001) as well as to the runaway growth limit of ~ 700 km 
derived from the proposed eccentricities and inclinations of 
the orbits of Fomalhaut's planetesimals (e.g., section 5). The 
distribution of runaway planetesimals is discussed in section 
7. 

We note that the shape of the primordial size distri- 
bution determines how the brightness of the disk evolves 
with time, and that this evolution can be constrained us- 
ing the distribution of the brightnesses of observed debris 
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Figure 8. Inferred size distribution of planetesimals in Fomal- 
haut's disk: A coUisional cascade extends from ~ 4 km planetes- 
imals down to 7 fim dust grains. Smaller dust grains are blown 
out of the system by radiation pressure. At sub-mm and shorter 
wavelengths, we only see those members smaller than 0.2 m. Plan- 
etesimals larger than ~ 4 km have an unknown primordial size 
distribution, however a reasonable distribution may be one that 
extends up to 1000 km with the same slope as the cascade (see 
section 5.3). Planetesimals larger than ~ 1000 km, if they exist, 
must have grown by runaway growth and would form a separate 
population (see section 7.2). 



disks. However we leave this discussion to a future paper. 
Our understanding of the size distribution of planetesimals 
in Fomalhaut's disk from this study is summarized in Fig. 0. 



6 INTRINSIC CLUMPINESS 

Collisions between planetesimals cause the disk to be in- 
trinsically clumpy, since after a collision the fragments of 
the target and the impactor move away from the impact 
site at a velocity determined by the kinetic energy imparted 
to them in the collision. This means that the debris from 
each planetesimal forms an expanding clump the centre of 
which follows the orbit of the parent planetesimal. At any 
one time we would expect to see clumps in the disk with 
a range of both fluxes (those from different magnitude col- 
lisions) and sizes (smaller ones that were created recently 
and older, more extended, clumps). The probability of our 
witnessing a coUisional clump of a given brightness and phys- 
ical size depends on: the rate of collisions between different 
sized planetesimals (section 6.1); the amount of dust such 
collisions produce, which depends on the size distribution of 
coUisional fragments (section 6.2); and on the rate at which 
this dust precesses around the orbit of the parent planetes- 
imal, which depends on the ejection velocity of coUisional 
fragments (section 6.3). 

6.1 Collision rates 

The number of coUisional events per unit time that occur in 
the disk between (target) planetesimals of size D to D+dD 
and those (impactors) of size between Dim and Dim+dDim is 
RcoiiD, Diin)dDiinn{D)dD, where the rate of collisions onto 
individual planetesimals, -Rcoi, was defined in equation (|ll|). 
These collision rates are most strongly affected by the disk's 
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size distribution and the total amount of cross-sectional area 
in it, but are largely unaffected by the eccentricities and in- 
clinations of disk material, as _Rcoi oc ^1 + 1.25e2/P. These 
orbital parameters do, however, have a strong effect on the 
outcome of these collisions. 



6.2 Collisional clump flux 

The size distribution of collisional fragments was not dis- 
cussed in BA99, since it was not permitted by the resolution 
of their models. Here we use the results of Campo Bagatin 
& Petit (2001; hereafter CPOl) which considers the size dis- 
tribution of collisional fragments expected from geometrical 
constraints. In the CPOl model, a planetesimal of mass M is 
broken up sequentially into smaller fragments, where the size 
of each fragment is determined by the shape of the volume 
remaining after the previously created fragments have been 
removed. The physical basis of this approach is that frag- 
mentation proceeds by the coalescence of flaws that propa- 
gate through the target after impact. The size distribution 
they flnd depends on the mass of the largest remnant cre- 
ated in the collision, fir = Mir/M. However, apart from the 
few largest fragments, the distribution was found to follow a 
power law with an index gc ~ 1.91 — 1.97, in agreement with 
their analytical model which predicted a distribution with 
gc = 1.93. This size distribution also concurs with experi- 
ments which flnd that while gc can fall anywhere in the range 
1.6-2.6 (Fujiwara et al. 1989; Giblin et al. 1998; Housen & 
Holsapple 1999), values close to 1.9-2 are more common. 
While it is not clear why geometrical arguments should ap- 
ply in the gravity dominated regime, since here each frag- 
ment is a gravitationally reaccumulated conglomerate of the 
smaller fragments created in the collision, CPOl proved that 
their model reproduces the size distributions of the aster- 
oid families, the members of which should be reaccumulated 
fragments (Campo Bagatin et al. 2001), and which have ob- 
served slopes of gc = 1.83-2.17. 

Here we assume that the size distribution of fragments 
smaller than the second largest remnant, that with diame- 
ter D2, follows a power law size distribution with gc = 1.93 
down to a lower limit -Dmin.ci. It may be inappropriate 
to extrapolate the CPOl size distribution to small frag- 
ments, since their distribution may be dominated by differ- 
ent physics to that of larger fragments. For example, Durda 
& Flynn (1999) found a knee in the size distribution of 
dust created by the break-up of composite materials due 
to flaws at material boundaries. This knee occurs at the size 
of crystals embedded in the planetesimal, (~ 1 mm in size 
in their experiment). A flattening of the size distribution for 
dust smaller than 1 mm was also reported by Fujiwara et 
al. (1989). In any case, we note that while fragments smaller 
than the radiation pressure blow-out limit, Dmin, may be 
created in collisions, they would be removed from the disk 
within an orbital timescale by radiation pressure. Thus we 

set -Dmin,cl ^ -Djxiin. 

Conserving the volume of the original planetesimal, and 
scaling the fragment distribution such that there are two 
fragments larger than D2 , we find that 



The total amount of cross-sectional area created by the 
break-up of a planetesimal is thus 

0-col(£',Am) = Cr„,i„CgJl -/lr(D, Am)]'''^^' 

x(D/A„i„,ci)'''="', (18) 

where C,, = [(2gc - 2)/(gc - 5/3)](2/gc - l)"'^-^ = 0.323. 
Equation ( [l8| ) is only valid when D2 < -Dir, i.e., when fi^ > 
1 — 0.5gc ~ 0.035. When /ir < 0.035 we assume that there 
is more than one fragment of size Dir, and that the gc = 
1.93 distribution applies for smaller fragments. Thus we get 
a new expression describing the outcome of such a super- 
catastrophic collision: 

Crcol(£',Am) = cr„,in[(2-gc)/(gc-5/3)]/i,(D,Am)''^"'' 
x(75/A.i„,ci)''^-'. (19) 

Note that the cross-sectional area given in equations ( [l8| ) 
and describe that of the debris produced by one of the 
planetesimals involved in the impact, that which we call the 
target, and that in both cases this cross-sectional area is 

min,cl 

BA99 found that the size of the largest remnant of 
a collision, whether rubble pile or intact fragment, is de- 
termined by the collisional energy in the following way: 
/ir = 0.5 — s((3/Qn — 1), where s is close to 0.5 for all materi- 
als and impact velocities, and Qo is the dispersal threshold 
at the relevant collisional velocity (i.e., not averaged over all 
impact parameters as before) . This is consistent with the re- 
sults of experiments of catastrophic impacts in the strength 
regime, which found that fi^ = 0.5(Qd/Q)^'^* (Fujiwara et 
al. 1989). It is also consistent with the results of experiments 
of cratering impacts in the strength regime: e.g., Holsapple 
(1993) used scaling arguments to show that in the strength 
regime, Mc.atcr/M = 1 - /i, cx (Q/t;c'oi)(«coi/Qs)'''^; for 
fj, — 2/3 and scaling to /h- — 0.5 when Q = Qg, we find that 
fir = l — 0.5Q/Qg. For collisions in which /ir > 0.5, we use 
the BA99 results with s = 0.5 so that 

Vcoi{D, All 



flr{D,Dlrn) = 1^0.5 



Dim 
Xtc{D)D 



vtc{D) 

■y(D) 



(20) 
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vtc{D) 


2. 48-1. 247(D) 


Am 




Vcol{D, Am) 





D2/D 



(1 - /ir) 



1/3 



(17) 



where we remind the reader that oc nj^j . For more en- 
ergetic collisions, however, we use the experimental results 
summarized in Fujiwara et al. (1989), since these cover a 
wider range of collisional outcomes, and BA99 only consid- 
ered collisions in which /ir ~ 0.5. Thus for /ir < 0.5, we 
use 

/lr(£',Am) =0.5 



which is consistent with the BA99 results in this regime. 
Here we assume that /ir averaged over all impact parameters 
can be obtained using the averaged QJ, model of section 5. 
We also set the physical constraint that /ir > (Dmin/15)^. 

Consider the dust clumps created in collisions between 
planetesimals in the Fomalhaut disk. In the following dis- 
cussion, the emission from these clumps is calculated from 
equation (^) assuming that the composition of this dust is 
the same as that inferred for the smooth ring from the SED 
modelling of section 4, and that this material is spread across 
the same range of distances from the star as the smooth 
ring (i.e., 150 ±25 AU). Fig. ^ shows the size of impactor 
required to produce debris clumps of target material that 
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10 10- 
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Figure 9. The size of impactor planetesimals in tlie Fomalhaut 
disk required to breals up tlie target into fragments that emit a 
given fiux, F^i, at 450 /im. The fragments are assumed to fol- 
low a power law with = 1.93 down to the radiation pressure 
blow-out limit of 7 fim. The solid and dashed lines correspond to 
the models for the coUisional properties of Ice and Basalt do- 
scribed in the text; in this target size range Weak Ice has the 
same properties as Ice. The four lines for each model correspond 
to -Fcl = 1 MJyi 10 iiJyi 0.1 nJy ^^<^ l pJy- The dotted lines de- 
mark the limits of threshold catastrophic collisions for the Ice 
and Basalt models. Collisions below these lines result in crater- 
ing from the target planetesimal, while those above result in the 
complete disintegration of the target. 



emit a range of clump fluxes at 450 ^m, assuming a frag- 
ment cut-off at -Dniin.ci ~ Dmin- This shows that for cra- 
tering events (i.e., for collisions below the dotted lines), the 
amount of dust produced depends mostly on the size of the 
impactor, not on the size of the target; i.e., whether an im- 
pactor falls on a planetesimal much larger than itself or one 
closer to its own size, a similar amount of dust is released 
from the target. Similarly, the amount of dust resulting from 
the catastrophic destruction of a target planetesimal (in a 
collision above the dotted lines) is largely independent of 
the size of planetesimal impacting the target. It is also in- 
teresting to note that debris from the smaller planetesimal 
in a cratering collision can be brighter than that from the 
larger planetesimal. 

If we had chosen the size distribution to be cut-off at 
I^min.ci ~ 1 mm, instead of 7 /im, then while the cross- 
sectional area resulting from a given collision would be de- 
creased by a factor of (1000/7)"'^ ™, the flux per unit cross- 
sectional area would have increased, because while the extra 
7 fim to 1 mm particles contain a lot of cross-sectional area, 
they don't emit efficiently at 450 fim. As the clump's 450 
fim flux in Jy is 6.5 x 10"'^ times its cross-sectional area in 
AU^ if the cut-off is at 7 /im, but is 0.11 times this area 
if the cut-off is at 1 mm, this means that with this larger 
cut-off, the fluxes in Fig. ^ should be reduced by a factor of 
~ 3. 

To determine what the clumps in Fig. ^, which are de- 
fined by their 450 fim emission, would emit at other wave- 
lengths, consider Fig. ^ which shows the SED of their emis- 
sion normalised to 30 mjy at 450 ^m (the flux of the ob- 
served clump). This SED is shown for different Dmin.ci as 
its shape is strongly dependent on this cut-off due to the 
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1000 



Figure 10. The SED of emission from a clump in Fomalhaut's 
disk assuming it is composed of grains with a size distribution 
with Qc = 1.93 that is cut-off at the small size end below -Dniin cl- 
The clump's SED is normalised to 30 mJy at 450 fim, which is the 
flux of the clump observed by Holland et al. (2002; see Fig. ^). 
Also shown is both the observed SED (symbols) and the fit to 
this SED assuming all this omission comes from the smooth ring 
(solid line). 



increased temperature of small grains (Fig. ^ . For example, 
collisions corresponding to 1 /xJy at 450 /xm for a clump with 
a cut-off at 7 /im (i.e., those close to the top right of Fig. ^ 
would result in emission of 11 /xjy at 25 fj,m and 43 fiJy at 60 
^m, while if fragments from the same collisions had a cut-off 
at -Dmin.ci = 1 mm, this would result in emission of 0.3 /j,Jy 
at 450 fim, just 7 njy at 25 fim, and 0.7 fiJy at 60 fim. Note 
that if the coUisional size distribution extends down to the 
blow-out limit, the clump would appear brighter relative to 
the smooth ring at mid- to far-IR wavelengths than in the 
sub- mm. 



6.3 Clump growth rates 

The rate at which the physical extent of a clump grows is 
determined by the ejection velocity of coUisional fragments. 
This is usually discussed in terms of the parameter /ke , the 
fraction of the impact energy, _Ecoi = 0-5M Mimv'^^i/ {M + 
Mim), that ends up as kinetic energy of the fragments, a 
parameter that is not well known. Laboratory experiments 
of impacts of cm-sized objects imply that the kinetic energy 
imparted to the largest fragments produced in the collision 
(i.e., those containing 70-80 per cent of the target mass) 
is a fraction /ke = 0.3-3 per cent of the impact energy 
(Fujiwara & Tsukamoto 1980), with an unknown quantity 
imparted to smaller fragments. Studies of the asteroid fam- 
ilies, on the other hand, imply that family-forming events 
had /ke ~ 0.1 (e.g., Davis et al. 1989). However, SPH sim- 
ulations of impacts in the gravity regime imply even higher 
values of /ke = 0.2-0.4 (Love & Ahrens 1996). It has also 
been reported that /ke increases as the velocity of the im- 
pact increases oc Q"'^ (e.g.j Fujiwara & Tsukatmoto 1980; 
Paolicchi et al. 1996), and it is expected that collisions with 
rubble pile planetesimals have a lower /ke than those with 
solid bodies due to the inefficient coupling of impact energy 
(e.g., Campo Bagatin et al. 2001). While clearly a simplifi- 
cation, in this paper we assume that /ke ~ 0.1, and that 
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this is true in both the strength and gravity regimes. We 
also assume that this energy is split evenly between the tar- 
get and the impactor, since this was found to be the case 
both for cratering events in the strength regime (Hartmann 
1988) and for collisions in the gravity regime (Love & Ahrens 
1996), although we recognize that the division of this energy 
depends on the impact parameter (Love & Ahrens 1996) 
among other things. 

The kinetic energy, /keScoI, is not evenly distributed 
across the fragments. Fragments are found to attain a range 
of velocities, the distribution of which can be approximated 
by the relation f{v) oc v~'', where f(v)dv is the fraction of 
material with velocities between i; and v+dv. The power law 
index, k, has been reported with values between 3.25 (Gault 
et al. 1963) and 1.5 (Love & Ahrens 1996), where values of 
k < 2 imply that most of the kinetic energy is carried away 
by a small mass of material which is moving very fast (e.g., 
Love & Ahrens 1996; BA99) and which usually originates 
near the impact site. There also appears to be a dependence 
of the velocity of a fragment on its mass in that, on average, 
lower mass fragments move faster than more massive ones. 
The mean velocity of fragments of size D deduced from lab- 
oratory experiments is reported to be oc , where r lies 
anywhere between (e.g, Giblin et al. 1998) and 0.5 (e.g., 
Nakamura & Fujiwara 1991), although the scatter about this 
mean tends to dominate over any strong trend (Giblin et 
al. 1998) . Very little of the kinetic energy is found to be im- 
parted to the largest remnant (Nakamura & Fujiwara 1991; 
Michel et al. 2001). 

Here we calculate the characteristic ejection velocity of 
target fragments assuming that the available kinetic energy 
is distributed among all fragments except the largest rem- 
nant; i.e., 0.5/ke£'co1 = 0.5(1 — fu-)Mv'^^, so that 

/KEQ(D,Am)/[l + (Am/D)'] 



by 



(22) 



[1 - MD, Am)] 

This means that target fragments created in a cratering 
event, wherein /ir is given by equati on (|2c|) a nd Dim ^ D, 
have an ejection velocity of Vej ~ y^/KEQ^j. Thus crater- 
ing debris from a given target has the same ejection ve- 
locity regardless of the size of the impactor. This is be- 
cause while the amount of kinetic energy imparted to frag- 
ments is larger for larger impactors, there is also a greater 
mass of material ejected from the largest remnant for this 
to be shared between. This relation also means that the 
ejection velocity is lowest for cratering events that occur at 
the transition between strength and gravity scaling, since 
these are the weakest planetesimals (Fig. ^). Equation ( |22[ ) 
also means that the ejection velocity of the fragments of the 
smaller planetesimal in the collision, for which /ir ~ 0, is 



fraction (22 per cent in this case) of its original velocity. 
This means that for planetesimals for which gravitational 
focussing is not important, all impactor debris has the same 
ejection velocity, which is ~ 90 m s~^ in the Fomalhaut disk. 

In the gravity regime, some of the kinetic energy im- 
parted to the fragments is converted into gravitational 
energy, since the fragments have to overcome the grav- 
ity of the largest remnant. The energy required to dis- 
perse all fragments except the largest remnant is -Egrav = 
1.2{GM^/D)[1 - /f/^]. Thus the characteristic velocity of 
the fragments once far from the largest remnant, Voo, is given 



2 



0.47rGpA[l-/i,(I>,Am) 
x[l-/i,(AA„0]'\ 



5/3, 



(23) 



(24) 



where this assumes that the kinetic energy given to the in- 
dividual fragments of the largest remnant is small. How- 
ever, using this model we find that for collisions between the 
largest planetesimals {D > 100 km) in the gravity scaling 
regime in Fomalhaut 's disk, Vcj is sometimes slightly smaller 
than This occurs when QJ, < {l/3)npGD^ / /ke- This 

discrepancy is due to our extrapolation of Qq from the BA99 
models to larger planetesimals and different impact veloci- 
ties. We already know that the collisional model may not 
be applicable to the break-up of planetesimals larger than 
~ 700 km diameter in this disk, and such collisions are dis- 
cussed in section 7. However, it may also be that /ke is 
higher in the gravity regime. Setting /ke = 0.34 would re- 
solve the discrepancy is this size range. A higher /ke may be 
expected because a much smaller fraction of the impact en- 
ergy is used in shattering the target; e.g., if /ke oc (Q/Qs)'^'^ 
(Fujiwara & Tsukamoto 1980), then threshold catastrophic 
collisions in the gravity regime would have a very high /ke 
(although some of this kinetic energy would go into the in- 
dividual fragments of the largest remnant before they reac- 
cumulate, at which point it would be converted into inter- 
nal energy). Here we account for this discrepancy by artifi- 
cially increasing /ke where required so that Vej > l.Oliigrav. 
The reasoning behind this is that in the gravity regime the 
fragments attain a range of velocities, and while those with 
Vej < Wgrav remain bound to the largest remnant, most of 
those that escape must have Voj just larger than Wgrav. 

The rate at which a clump appears to grow depends on 
the rate at which the fragments precess around the orbit of 
the parent planetesimal. The radial and vertical structure 
of the clump change with time as well, however the spatial 
extent of these variations are at or below the level of the az- 
imuthal growth. A fragment's precession rate is determined 
by the velocity with which it is ejected in the direction of 
the orbital motion, dvk, since this causes a change in the 
semimajor axis of the fragment's orbit, and hence in its or- 
bital period, diper/ipor = —3dv]^/v]^. The ejection velocity 
field is generally found to be isotropic both in experiments 
(Giblin et al. 1998) and analysis of the orbital elements of 
the asteroid families (Zappala et al. 1996) . Thus the average 
velocity of all fragments in the direction of orbital motion is 
d^k ~ Voo/V^, and the rate of growth of the total azimuthal 
extent of the clump is 



\/0.5/KEt'coi(-D, Am); i.e., the impactor retains a fixed g — 2\/3uoo/ffl. 



(25) 



Before continuing we define the characteristic lifetime 
of a clump of target fragments created in a collision be- 
tween the target (of size D) and an impactor (of size Am), 
tci{D, -Dim), to be equal to the average time it takes for the 
clump fragments to occupy half the ring (i.e., Tr/9): 



fcl(D, Am)/V.- = (l/4%/3)KAoo(D, Am)]. 



(26) 



Thus in this model impactor debris from collisions between 
planetesimals in Fomalhaut's disk remains in a clump for 
about 6 orbital periods, while target debris clumps last for 
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between 10 and 700 orbital periods, depending on the tar- 
get's dispersal threshold (clumps from weak, 10-1000 m, 
planetesimals last longest). 

6.4 Discussion 

Combining the collision rate from section 6.1 and the clump 
growth rate from section 6.3, we find that the number of 
clumps we expect to see at any one time that are smaller in 
azimuthal extent than 9, and that were created in collisions 
between (target) planetesimals of size D to D+dD and those 
(impactors) of size between Dim and Dim+dDim is given by 

iV(< e) = iicoi(D, Din.)dDimri{D)dD{e/e). (27) 

We also know how bright those clumps should be, since this 
was defined in section 6.2 (Fig. Thus we can integrate 
over all impacts that create clumps brighter than a specific 
limit to get the total number of clumps we expect to see 
in the disk that are smaller in azimuthal extent than 9 and 
brighter than Fci, A*'(< 6^, > -Fci)- 

Clumps that are greater than ~ 180° in extent are no 
longer clumps, but radially confined rings that blend into the 
background disk. Furthermore, the emission from such rings 
could be affected by the subsequent coUisional evolution of 
their constituents^ Here we define the clump function to 
be the number of clumps we expect to see in the disk that 
are smaller than 180° in extent and brighter than Fd, N{9 < 
180°, > -Fci)- The number of clumps brighter than this limit, 
but with a maximum extent smaller than 180° , can then be 
determined using the relation: 

N{< 9, > Fci) = N{< 180°, > Fci) X (9/180°). (28) 

In Fig. |l^ we have plotted the clump function in Fo- 
malhaut's disk assuming that the disk is comprised of: (i) 
just the planetesimals in the coUisional cascade; (ii) the col- 
lisional cascade plus an additional ~ IMj component of pri- 
mordial planetesimals with a size distribution defined by 
Qp = 1.84 extending up to those 1000 km in diameter. The 
steepness of the clump function in this plot, combined with 
the linear growth rate of clumps, implies that this clump 
function can be loosely interpreted as the number of clumps 
we expect to see of fiux ~ F^i that are 90-180° in extent. 
The origin of these clumps could be in any collision on the 
appropriate line in Fig. ^ However, we find that clumps 
produced in the gravity regime are most likely to be target 
debrisj^ from a threshold catastrophic collision; e.g., any 1 
nJy clump we observe is most likely material created in a 
collision between two bodies ~ 300 km in diameter. 

We define the intrinsic dumpiness of a disk as the 
flux at which the clump function, N{< 180°, > Fd) = 1. 
This can be interpreted as the approximate fiux of the 
largest single clump we would expect to see in the disk. 
However, as a disk contains a distribution of clump sizes, 

* The subsequent coUisional evolution of coUisional fragments 
would not affect their emission by much until 6 > 180°, as the 
coUisional lifetime of this material, assuming it is the same as that 
of the background disk (Fig. is much longer than the age of 
the clump (equation This evolution would tend to reduce the 
clump's brightness, as collisions remove the smallest dust first. 
^ As target and impactor debris remains in a clump for different 
lengths of time, we have treated the two components separately. 




Clump flux at 450 f^rn, F^,, Jy 



Figure 11. The clump function, N{e < 180°, > Fcl), of the Fo- 
malhaut disk, i.e., the number of clumps we expect to see that are 
brighter than but smaller than 180° in azimuthal extent. The 
solid, dotted and dashed lines correspond to the models for the 
coUisional properties of Ice, Weak Ice and Basalt described in 
the text. The two lines for each model assume that the size distri- 
bution with = 1.84 extends up to maximum planetesimal sizes 
of 4 km and 1000 km. The intrinsic dumpiness of the disk (the 
Fci at which N{0 < 180°, > Fd) = 1) is shown by the horizontal 
long dashed line. Its observed dumpiness, F^ = 30 mjy, is shown 
by the vertical long dashed line. For collisions to be the likely 
cause of this dumpiness we would expect the clump function to 
pass through the intersection of these lines. 



smaller fainter clumps would also be observed down to a 
level determined by the beam size. Observations with a 
beam of width FWHM (in arcsec) , would see all clumps with 
9 < 0PWHM, where SpwHivi = 360° x (FWHM r^,)/{2nrd) 
(rci is the distance of the clump from the star in AU) as un- 
resolved point sources; those with Spwhm < 9 < 180° would 
be resolved, while larger ones would be rings rather than 
clumps. Thus, the number of clumps we would expect to see 
as point sources above a certain flux would be Sfwhm/180° 
times the clump function, and the brightest unresolved point 
source we could expect to see would be the flux at which the 
clump function, N{< 180°, > Fd) = 180°/6Ifwhm. 

Fig. ^ shows that if Fomalhaut's disk is comprised 
solely of planetesimals in the coUisional cascade, its largest 
(4 km) planetesimals are so numerous, and their clumps last 
so long, that at any one time we would expect to see roughly 
one million clumps resulting from their destruction that are 
< 180° in extent. Thus even if the sensitivity of our obser- 
vations was sufficient to detect these clumps (which would 
emit just ~ 10 pJy), the disk would appear smooth unless 
our beam size was also < 60 ^arcsec! However, if the disk 
also contains a primordial planetesimal distribution as de- 
scribed in (ii) above, the disk's intrinsic dumpiness would be 
~ 2 /ijy (Fig. ^l|). Brighter clumps could also exist in the 
disk (up to 10 /iJy), however such large coUisional events 
happen infrequently. For these clumps, the clump function 
defines the probability of our witnessing a clump that bright 
and smaller than 180°; i.e., at any one time we have a 1:100 
chance of seeing a 10 nJy clump (if the proposed distribution 
is correct). 
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Could intrinsic dumpiness be the likely cause of the 
clump observed in Fomalhaut's disk (section 2)? If so we 
would expect the derived intrinsic dumpiness at 450 
to be 30 mjy; i.e., the clump function should pass through 
N{< 180°, > 30 mJy) « 1. This is not true in either of the 
examples shown in Fig. ^ However, we have no reason to 
expect that qp = 1.84, nor that this distribution extends 
up to -Dmax = 1000 km. A lower would increase the num- 
ber of 1000 km planetesimals in the disk, and so increase the 
number of clumps expected from their destruction (although 
this would also imply a disk more massive than IMj, equa- 
tion (|l^)). The clump function could also be underestimated 
due to uncertainties in our model. Particularly, the lifetime 
of clumps produced by the destruction of planetesimals in 
the gravity regime is not well known (as discussed in sec- 
tion 6.3), and may have been underestimated. These factors 
could mean that the intrinsic dumpiness is actually at the 
level of the break-up of the largest planetesimal in the disk. 
However, if Dmax ~ 1000 km, this is just ~ 10 /iJy, some 
3000 times lower than that required for intrinsic dumpiness 
to explain the observed clump. 

However, it is possible that our model also underes- 
timates the amount of dust produced in collisions in the 
gravity regime, since smaller dust grains might be predom- 
inant in the escaping fragments from these collisions. If the 
velocity imparted to the fragments is size dependent, and 
smaller fragments attain higher velocities, then it would be 
the largest fragments of the original distribution that reac- 
cumulate, while the smaller ones are more likely to be re- 
leased. Furthermore, the surfaces of large planetesimals are 
expected to be covered with large quantities of small dust 
grains in the form of a regolith (see section 7.2); these grains 
would be preferentially released in a collision. It has not 
been possible to include such effects in our model because 
the outcome for /im-sized dust from km-scale impacts has 
yet to be studied. Based on arguments outlined in section 
7.1, the most our model could be underestimating the flux 
by is a factor of 1000^. However, the real value is likely to 
be much less than this. In any case, collisions among non- 
gravitationally focussing planetesimals result in insufficient 
flux for them to be the cause of the observed clump. 

6.5 For Future Consideration 

Could Fomalhaut's intrinsic dumpiness be observed with fu- 
ture instrumentation? With a beam size of 0.6 arcsec. Fig. |l^ 
implies that we could expect to observe clumps in Fomal- 
haut's disk between ~ 0.1 /iJy (unresolved point sources) 
and 2 /xjy (resolved clumps) at 450 fim. Depending on the 
amount of /xm-sized dust produced in collisions, these lim- 
its correspond to dumps of between 1 /xjy and 0.02 mJy 
at 25 /im (Fig. |l^. A 0.6 arcsec beam is close to that of 
the ALMA in compact configuration at 450 fim, or that of 
a mid-IR instrument operating at 25 /xm on an 8 m tele- 
scope. While the dumpiness of Fomalhaut's disk proposed 

^ The collisional model predicts that the destruction of a 1000 
km planetesimal would result in dust that emits a flux of 10 ^Jy, 
whereas the maximum possible emission from a pulverised 1000 
km body, i.e., when it is ground into 150 fim—l mm dust, is 10 
mJy (Fig. ^) 



above is at least an order of magnitude below the anticipated 
sensitivities of such instruments, the above figures include 
uncertainties (such as the potential factor of 1000 underes- 
timate of the amount of dust created in a collision, section 
6.4) which could render these clumps detectable. If so, their 
observed magnitude and size distribution could be used to 
constrain the primordial planetesimal population. 

Also, we note that this type of dumpiness could provide 
confusion for a mission to find terrestrial planets using space- 
based 10 fim interferometry, such as the designs proposed 
for DARWIN or TPF (e.g., Mennesson & Mariotti 1997). 
The biggest potential impact of disk dumpiness on such a 
mission would come from a warm zodiacal cloud-like disk ^ 
3 AU from the star; this topic is discussed in Wyatt (2001). 
However, here we point out that if viewed at the distance 
of Fomalhaut, the Earth would have a 25 /im flux of ~ 0.8 
/ijy, while the primordial planetesimal distribution proposed 
in Fig. |l^ implies that the Fomalhaut disk could contain a 
comparable ~ 0.2 /ijy clump that would be unresolved in a 
0.05 arcsec beam, which would be that of a 100 m baseline 
interferometer at this wavelength. 



7 THE COLLISIONAL CLUMP HYPOTHESIS 

Consider the magnitude of collision required to produce the 
clump observed in the Fomalhaut disk. Section 6.4 showed 
that this clump cannot have been produced by the intrinsic 
dumpiness of the collisional cascade, or by the break-up of 
the primordial planetesimal population, because such colli- 
sions do not produce enough dust. Thus we start by deter- 
mining the smallest planetesimal that could have produced 
the clump using a simple model in which the clump is com- 
prised of grains of just one size. 

7.1 Single size dust grained clump 

If the dust grains in the clump are all of the same size, Dd, 
then Fig. |l2p , shows how the observed 450 fim flux of the 
clump, F^(450 /im) = 30 mJy, determines the total cross- 
sectional area of material in the clump, ad (found by in- 
verting equation for different D^i- Thus the projected 
area of this clump must be at least 0.2 AU'^, much larger if a 
signiflcant fraction of the material is in dust smaller than 1 
mm. Fig. ^b also shows how this flux determines the total 
mass of material in the dump, where md — (2/3)p(Tci-Dci- 
Thus the clump must have a mass of at least a thousandth 
that of the Earth, and would only be this small if all the 
grains in the clump are between 150 fj,m and 1 mm in di- 
ameter. If there is a signiflcant fraction of mass in particles 
smaller than 150 ^im then we would infer a more massive 
clump, since these particles do not emit efficiently at 450 
^m. Similarly we would infer a more massive clump if there 
is a significant fraction of mass in particles larger than 1 mm. 
It would be possible to constrain the size of particles in the 
clump if the clump could be imaged at shorter wavelengths 
(see Fig. ^c). However, for now we can rule out the possi- 
bility that the clump consists of grains that are all smaller 
than ~ 100 /xm, since the emission from such a clump would 
exceed that observed, presumably from the smooth ring, at 
mid- to far-IR wavelengths. 

Thus if the material in this clump originated from a 
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Figure 12. The total cross-sectional area (a) and mass (b) of 
material in Fomalhaut's clump, as well as this clump's SED (c), 
assuming it is comprised of grains of just one size, D^i, and setting 
the clump's emission to 30 mjy at 450 /im. As in section 6, the 
clump is assumed to be composed of the same material as that 
inferred for the smooth ring and the clump's emission is calculated 
assuming that this material is spread across the same range of 
distances from the star as the smooth ring. Also shown in (b) is 
the diameter of the body that the clump mass would correspond 
to when recombined into a single body, and (c) also shows the 
observed SED (symbols) and the fit to this SED assuming all this 
emission comes from the smooth ring (solid line). 



single body, that planetesimal would have to have been at 
least 1400 km in diameter. To understand how such a large 
body, which must have formed through runaway growth, 
and for which gravitational focussing is important, interacts 
with the planetesimal disk, we start by considering how such 
a planetesimal grew in the first place. 



7.2 Runaway planetesimal growth 

The runaway population are the largest members of the dis- 
tribution of planetesimals growing in the protostellar nebula. 
They grew initially in the same way as the rest of the plan- 
etesimals, by the accretion of nearby, usually much smaller 
planetesimals. However, kinetic energy exchange in encoun- 
ters with other planetesimals (i.e., dynamical friction) led to 
larger planetesimals having more circular orbits (Wetherill 
& Stewart 1993). The consequent low relative velocity and 
high gravitational focussing factor for collisions with large 
bodies (e.g., equation (|l^)), meant that large bodies (once 
larger than ~ 1 km) grew more rapidly than smaller plan- 
etesimals, creating a bimodal size distribution - an accumu- 
lation cascade and a population of large, so-called runaway, 
planetesimals. This runaway growth would have been rela- 
tively short, typically several Myr (e.g., Kenyon 2002), and 
would have been halted once the runaways became large 
enough to perturb the disk just enough to affect their coUi- 
sional velocities (i.e., once larger than a few hundred km). 
The runaways then continue to grow, but at a slower rate 
than smaller planetesimals, so-called oligarchic growth (e.g., 
Kokubo & Ida 2000). Once these runaways are large enough, 
however, their gravitational perturbations increase the im- 
pact velocities of members of the accumulation cascade so 
that collisions result in their net destruction, rather than in 
their net accretion. 

During this oligarchic growth, the runaways try to re- 
main dynamically isolated from each other. However, this is 
only possible as long as their separation is at least 5 Hill's 
radii (Kokubo & Ida 1995), where a planetesimal's Hill's ra- 
dius is r-H = a{Mpi/2>Mi,y'^^ . A runaway's increasing mass 
means that any isolation is not permanent and interactions 
with other runaways are inevitable. These interactions would 
result in both collisions amongst the runaway population, 
as well as increases in the eccentricities of this population, 
thus increasing their feeding zones (e.g., Weidenschilling et 
al. 1997) and making mutual collisions more likely. Their 
growth would eventually be halted once they had scattered 
or accreted all of the planetesimals in their feeding zone. 
The angular momentum exchange caused by this scattering 
could cause the orbit of the newly formed planet to migrate 
either towards or away from the star (e.g., Hahn & Malhotra 
1999). 

A runaway body would have started as an incoherent 
body built up of layers of regolith. Mutual collisions between 
runaways would have resulted in the runaway often being 
completely fragmented, then reaccumulating. Since smaller 
fragments attain higher ejection velocities, the reaccumu- 
lated body would be sorted according to grain size with the 
larger grains towards the centre and fine grains on the sur- 
face (e.g., Britt & Consolmagno 2001). Impacts of smaller 
bodies onto the runaway would result in the comminution 
of the runaway's surface regolith layer (e.g., McKay, Swin- 
dle & Greenberg 1989). Planetesimals in the solar system 
have such regoliths: the Moon's regolith is 5-10 m deep in 
the mare region, and is comprised mostly (80-90 %) of dust 
grains smaller than 1 mm in size (e.g., Heiken 1975), while 
that of 20 km diameter Phobos is thought to be more than 
100 m thick in places (Veverka & Thomas 1979), and the re- 
golith layer on 200 km-sized asteroids is predicted to extend 
to a few km in depth (Housen & Wilkening 1982). Thus we 
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expect runaways to be composed of a megaregolith structure 
on top of which is a deep dusty regolith layer. 

7.3 Observing the final throes of planet 
formation? 

Within this scenario there are three ways in which the ob- 
served dust clump could be associated with a runaway plan- 
etesimal in the ring: (i) it is unbound material generated 
in a recent collision between two runaways; (ii) it is bound 
material constantly replenished by impacts into a growing 
runaway; (iii) it is an infalling envelope of dust that is be- 
ing accreted from the collisional cascade onto the runaway. 
The clump could also be associated with a planet orbiting 
interior to the ring: (iv) its dust could have been created by 
the break-up of planetesimals that are trapped in resonance 
with the planet. 

7.3.1 Unbound runaway-runaway collisional clump 

If we assume that the M0/IOOO of dust that we see origi- 
nated in just one collision, then it must have originated from 
a parent planetesimal more than 1400 km in diameter. For 
this dust to be released from the gravity of the parent body, 
the impacting planetesimal must have had enough mass that 
the kinetic energy imparted to the fragments exceeded the 
escape velocity. Since impacts with gravitationally focussing 
planetesimals take place at the escape velocity, at best an 
impactor can only release /ke times its own mass. Thus the 
observed clump, if the product of one collision, must have 
involved the collision of two runaway planetesimals at least 
1400 km in diameter. In fact, the planetesimals may have to 
have been significantly larger than this, since experiments 
show that transmission of kinetic energy is much less efficient 
for collisions into sand than into rock (e.g., Holsapple 1993) 
leading some authors to speculate that /ke ~ 0.01 for colli- 
sions with such rubble piles (Campo Bagatin et al. 2001). 

How likely are we to be observing a collision between 
two runaways? Planet formation simulations do predict the 
formation of a significant population of runaways. Also, we 
know that runaway-runaway collisions do occur: based on 
the similarity of the abundances of isotopes of oxygen in 
lunar samples to those on the Earth, Wiechert et al. (2001) 
concluded that the Mars-sized impactor that collided with 
the Earth to form the Moon formed at the same heliocentric 
distance as the Earth (i.e., it was a nearby runaway). To 
have grown to > 1400 km, the runaways in Fomalhaut's 
disk must have grown in isolation from each other, with at 
least 5rH separation. This means that there can be at most 
one hundred M^/IOOO runaways in Fomalhaut's 50 AU wide 
ring. Since each collision would result in the merger of the 
runaways and so in the loss of one of the runaways, this 
leaves at most about 100 possible collisions resulting in the 
observed clump, that is if the runaways are not scattered 
out of the disk by their mutual perturbations before they 
collide. 

The length of time this material would remain in a 
clump depends on its ejection velocity. The only studies 
that can throw light on the outcome of such collisions are 
the SPH simulations of the formation of the Moon (e.g., 
Cameron 1997; Canup & Asphaug 2001) and of the aster- 
oid families (Michel et al. 2001). However, even these do not 



predict the outcome for collision fragments that are smaller 
than their resolution (> km), or indeed what would happen 
to the dusty regoliths on the colliding bodies. What they do 
show, however, is that in a collision between two runaways 
a significant mass of fragments does escape their gravita- 
tional field (see table 1 of Canup & Asphaug 2001). Thus 
such collisions would result in the formation of an unbound 
clump. Also the ejection velocities of fragments from col- 
lisions between massive, but non-gravitationally focussing, 
planetesimals imply a clump lifetime (equation (^^) of 20- 
30 orbital periods for collisions in the asteroid belt (Michel 
et al. 2001). Thus if this clump lifetime is also typical for 
runaway-runaway collisions in Fomalhaut's disk, we could 
only expect to see a clump from their mutual collisions for 
a total of ~ 3 Myr of the disk's life. 

This is an extremely short window to be witnessing the 
product of runaway-runaway collisions given the age of the 
system (3 Myr is 1-2 per cent of the disk's life). However 
the lifetime of a clump created in a collision between gravi- 
tationally focussing bodies in Fomalhaut's disk could be an 
order of magnitude longer than those of clumps created in 
the asteroid belt. Also planet formation models predict that 
the growth of the runaways to sizes large enough both to 
cause the observed clump and to ignite of the collisional 
cascade can take a few hundred Myr at the distance of Fo- 
malhaut's disk (e.g., Kenyon & Bromley 2001). Thus we 
leave this topic to be explored using more detailed planet 
formation and collision models, but remark that it appears 
very unlikely at this stage that runaway-runaway collisions 
are the cause of the observed clump. 

7.3.2 Bound collisionally replenished envelope (CRE) 

Impacts into a runaway planetesimal would launch a cloud 
of regolith dust from its surface. The majority of this dust 
would then collapse back onto the runaway. The maximum 
distance out to which dust would remain bound to the run- 
away would be to its Hill's radius, since outside this its grav- 
ity would be weaker than that of the star; this corresponds to 
~ 1.2 AU for a IM^ planetesimal orbiting at 150 AU around 
Fomalhaut. Thus for the observed dust clump, which has a 
projected area of > 0.2 AU'^ (section 7.1), to be bound to its 
parent planetesimal, that planetesimal would have to have 
a mass of at least O.OIM®. Here we consider whether Foma- 
lhaut's clump could be caused by repeated collisions onto 
a runaway from non-gravitationally focussing planetesimals 
sustaining a dusty envelope of regolith material extending 
to the Hill's radius. Since all collisions with runaways occur 
at the escape velocity, and the ejection velocity required for 
dust to reach to the Hill's radius is close to (just 5 x 10~^ 
times less than) that escape velocity, energy constraints im- 
ply that the maximum amount of regolith dust that an im- 
pacting planetesimal could launch to such a distance would 
be equal to its mass. Thus the maximum mass of material 
in such an envelope at any one time is equal to the mass 
accretion rate onto the runaway times the maximum length 
of time that the bound dust could remain in the envelope. 

The mass accretion rate of a runaway of size Drun is 
given by 

Maccr(Aun) = ( (7rp/6)ALi?col(Aun, Am)dAm. (29) 
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planets could not explain the observed dust clump. However, 
the flux of these CREs could surpass that of the runaway 
itself by up to a factor of 10,000 for the 450 flux of a 10 
planet (or perhaps even higher at shorter wavelengths, 
see Fig. p^) . These CREs could also be bright enough that 
it is the distribution of massive runaways that determines 
a disk's intrinsic dumpiness, rather than the distribution of 
non-gravitationally focussing planetesimals (section 6.4). 

7.3.3 Inf ailing dust envelope 

Another reason why there would be an enhanced density of 
dust in the vicinity of a runaway is that it is accreting dust 
from the disk's coUisional cascade. The total cross-sectional 
area of material in this infalling envelope is approximately 
the runaway's rate of accretion of cross-sectional area, given 
by 

aaccr(Drun) = /" (7r/4) A^m^col ( Aun , Am )d Am , (30) 

multiplied by the infall timescale, which is the freefall 
timescale from its Hill's radius (one quarter of an orbital 
period). The flux from this infalling material can then be cal- 
culated by considering that the disk emission of 0.595 Jy at 
450 iJ,m (Table 1) is predicted to originate from crtot = 33.7 
AU^ of material (section 4.3). The resulting flux is plotted 
in Fig. |l^, which shows that such an infalling envelope can 
be ignored as it is fainter even than the runaway. 



Figure 13. The interaction of a runaway planetesimal with the 
non-runaway population in Fomalhaut's disk: (a) the rate of ac- 
cretion of mass onto the runaway; (b) the flux of clump expected 
from this interaction due to both the coUisionally replenished en- 
velope (CRE) of regolith material and the infalling dust envelope, 
as well as the emission from the runaway itself (assuming a tem- 
perature of ~ 45 K). Since a runaway accretes most of its mass 
from the largest members of the population, the mass accretion 
rate and the flux from the CRE are shown for both Z?max = 4 km 
and 1000 km (dotted and solid lines respectively). 



This is plotted in Fig. |l3| a for size distributions that ex- 
tend up to 4 km (i.e., the coUisional cascade) and 1000 km 
(i.e., including a hypothetical primordial planetesimal pop- 
ulation) . Note that since the growth rate of runaways in the 
current disk is relatively slow, they must have formed when 
the eccentricities and inclinations of disk material were much 
lower, since then the accretion rate would have been higher 
(Aifaccr oc [//(e, /)] ~ ^ ) . The maximum amount of time dust 
from any one collision would remain in the envelope would 
be twice the free-fall timescale from the Hill's radius (e.g.. 
Love & Ahrens 1996), i.e., about half an orbital period, or 
650 years at 150 AU from Fomalhaut. Some material could, 
however, remain for longer if injected into orbit around the 
runaway (e.g., Canup & Asphaug 2001; Michel et al. 2001), 
The maximum 450 pim flux a given mass of material in 
Fomalhaut's disk could emit can be calculated from Fig. p^b. 
If all this mass was in dust 150 pim to 1 mm in size, this 
would emit 30 mJy/O.OOlM®. Fig. [l^ uses this to show the 
maximum flux of such a coUisionally replenished envelope 
(CRE). This shows that even the CREs surrounding lOMe 



7.3.4 Dusty resonant ring 

While no planets have been detected around Fomalhaut, the 
existence of a planetary system orbiting within the observed 
dust ring would naturally explain the lack of dust emission 
there, since such a system would have cleared the plan- 
etesimal disk within which it formed. As stated in section 
7.2, this clearing could have caused the planets to migrate 
outwards. In the course of this migration, the outer mean 
motion resonances of the outermost planet of the system 
would have swept through the disk of planetesimals out- 
side the planet's orbit. Many of these planetesimals would 
have become trapped in these resonances and would then 
have migrated out with the planet, while those that did not 
get trapped would have been scattered by the approaching 
planet. This is the mechanism which is invoked to explain 
the large number of Plutinos, Kuiper belt objects that are 
trapped in Neptune's 2:3 resonance (e.g., Jewitt 1999). Hahn 
& Malhotra (1999) estimated that the required migration 
of Neptune's orbit from 23 AU to 30 AU could have been 
caused by the clearing of a 50 M9 planetesimal disk over 
the first 50 Myr of the solar system. We note that these last 
two figures are comparable to the inferred mass and age of 
Fomalhaut's disk, thus it could be argued that similar mi- 
gration and trapping are to be expected in this system had a 
planet formed in a (now scattered) planetesimal disk interior 
to the one we see today. 

The paths of eccentric resonant orbits when plotted 
in a frame co-rotating with the planet's mean motion ex- 
hibit loops which occur at the pericentres of the orbits (see 
fig. 8.4c of Murray & Dermott 1999). As only those resonant 
orbits that have their conjunction with the planet when at 
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their apocentres are stable (at least when the planet has 
a low eccentricity), these loops only occur at specific loca- 
tions relative to the planet. Planetesimals that are evenly 
distributed within these stable resonances would be most 
densely concentrated in these loop regions. Thus, assum- 
ing that the fragments created by the collisional break-up of 
such resonant planetesimals also remain in resonance, the re- 
sulting dust ring would be clumpy (Wyatt 1999). The num- 
ber and distribution of these clumps would be determined by 
the fraction of planetesimals that arc trapped in the differ- 
ent resonances. The one clump in Fomalhaut's ring implies 
that the material is trapped in a 1:2 resonance, although 
other clumps could be hidden in the lobes. Very few Kuiper 
belt objects are observed to be trapped in Neptune's 1:2 
resonance, however this may be an observational limit, or 
because the primordial Kuiper belt did not extend out that 
far, and/or Neptune's orbit did not migrate far enough. The 
fraction of planetesimals that we require to be trapped in 
the the 1:2 resonance to cause the observed 30 mjy clump 
would be of the order of 5 per cent, since resonant material 
can spend a significant fraction of its time in the clump re- 
gion (the actual fraction is determined by its eccentricity) 
and while in the clump this material would be closer to the 
star and so hotter and brighter than that in the rest of the 
ring. 

If this interpretation is correct we would predict that the 
planet orbits at 80 AU with a period of about 500 years, and 
that it is currently located on the opposite side of the star 
from the dust clump. We would also predict that the clump 
would orbit the stax with the planet (direction unknown) 
causing motion on the sky of up to 0.2 axcsec/year, and 
that the clump would be hotter than the rest of the ring. A 
model for this resonant structure, and the constraints it sets 
upon the planet's mass, and that of the original planetesimal 
disk, is beyond the scope of this paper and will be discussed 
in a forthcoming paper. 

We note that the resonant ring we propose here is 
formed in a different manner to that proposed to explain 
the Vega and e Eridani dust clumps (Ozernoy et al. 2000; 
Wilner et al. 2002). The dust in their rings is trapped into 
planetary resonances through inward dust migration caused 
by P-R drag, rather than by outward planetary migration. 
Thus their rings are akin to the Earth's resonant ring (Der- 
mott et al. 1994). We contend that this mechanism would 
be inefiicient in this system, as grains would be destroyed 
by collisions before they migrate into the resonances by P-R 
drag (e.g., section 3.3; WDT99; Wyatt 1999). 



8 CONCLUSIONS 

We have presented a detailed study of collisional processes 
across the whole size range in material in the Fomalhaut 
debris disk. This is the first such study of a debris disk that 
we arc aware of. We have shown that the spatial distribu- 
tion of the dust derived from the 450 /im image of the disk, 
combined with modelling of its SED, are sufiicient to con- 
strain the dust's size distribution but not to set stringent 
constraints on its composition. The size distribution we in- 
ferred is the same as that expected from a disk in collisional 
equilibrium, although we recognize that our parameterizar 
tion of this distribution is only an approximation. The emis- 



sion that we observe originates from material in the size 
range 7 /im-0.2 m, however collisional lifetime arguments 
imply that this material originated in a cascade involving the 
break-up of planetesimals up to 4 km in size. The inferred 
mass of material in this collisional cascade is 20 — 30Af®. 
We speculate that there may be a population of primordial 
planetesimals with sizes between 4-1000 km, however these 
collide too infrequently for their collisional debris to have 
contributed to the dust we observe today. The slope of the 
size distribution of these planetesimals could be similar to 
that of the cascade and extend up to 1000 km based on mod- 
els of planet formation. Dust clumps caused by the break-up 
of these planetesimals would be at too low a level to be de- 
tected with current technology, however it may be possible 
to observe these clumps in the future (e.g., with the ALMA). 
With such observations we would be able to constrain the 
population of these large planetesimals — indeed this may 
be the only way of detecting their presence. There may also 
be a population of runaway planetesimals larger than 1000 
km which would still be growing. It may be possible to de- 
tect these directly, since they would always be enveloped by 
a cloud of dust launched from their regolith by repeated col- 
lisions. These dust envelopes could also be at a level which 
it may be possible to observe in the future. 

We propose two possible origins for the 30 mJy clump 
that is observed in Fomalhaut's disk. First it could have been 
created in a collision between two runaway planetesimals. 
We are limited in being able to assess the likelihood of this 
by uncertainties in the outcome of collisions between such 
massive bodies, as well as by uncertainties in the number 
and collision frequencies of such bodies at the end of planet 
formation. Our best estimate, however, shows that this pos- 
sibility is unlikely unless both the formation of the runaways 
and the ignition of the collision cascade occurred within the 
last few Myr. Another possibility is that ~ 5 per cent of the 
planetesimals in the ring were trapped in 1:2 resonance with 
an inner planet when it migrated out due to the clearing of a 
residual planetesimal disk. This hypothesis is strengthened 
by the fact that the mass and age of Fomalhaut's disk are 
comparable to that required for Neptune's migration and the 
consequent resonant trapping of Kuiper belt objects. Con- 
straints on the mass of Fomalhaut's hypothetical planet will 
be addressed in a future paper. We predict that the orbital 
motion of such a resonant clump (0.2 arcsec/yeax) could be 
detectable within a few years, and that this clump would be 
more prominent at shorter wavelengths. 

It is clear from this study that collisions are fundamen- 
tal to the formation and evolution of debris disks, not only as 
the mechanism which creates the dust we see, but as a cause 
of potentially observable structure in the disk. Unless future 
collisional and planet formation models can prove otherwise, 
it appears unlikely that Fomalhaut's observed clump is colli- 
sional in origin. However, this does not rule out a collisional 
origin for the clumps seen in other debris disks. These disks 
will be the subject of a future paper. 
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